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Projection of Bloch states obtained from quantum-mechanical calculations onto atomic orbitals is 
the fastest scheme to construct ab-initio tight-binding Hamiltonian matrices. However, the presence 
of spurious states and unphysical hybridizations of the tight-binding eigenstates has hindered the 
applicability of this construction. Here we demonstrate that those spurious effects are due to the 
inclusion of Bloch states with low projectability. The mechanism for the formation of those effects 
is derived analytically. We present an improved scheme for the removal of the spurious states 
which results in an efficient scheme for the construction of highly accurate ab-initio tight-binding 
Hamiltonians. 


I. INTRODUCTION 

The tight-binding method, even in its simplest imple¬ 
mentation, is a useful tool in the study of the electronic 
structure of molecules and solids.l^ The advantage of the 
method is the tractable and intuitive understanding it af¬ 
fords by distilling the electronic structure of complex sys¬ 
tems into physically transparent Hamiltonian matrices ex¬ 
pressed on a minimal basis set of atomic orbitals (AO). For 
realistic materials, the tight-binding (TB) matrix elements 
have been typically calculated by htting to experiments or 
higher levels of theory. The resulting models were success¬ 
fully applied in a panoply of complex materials with large 
supercell^^M^ and for problems where a localized basis sets 
are essential.l^ A major shortcoming of the tight-binding 
approach is the demanding fitting procedure that limits 
the application of the approach to well known materials 
and hinders the transferability of the parameters to bond¬ 
ing environments outside the assumed training set. 

In recent years, the accuracy and reliability of the TB 
models have been largely improved with the introduction 
of ab-initio tight-binding Hamiltonians. Here, the Hamil¬ 
tonian resulting from a fully self-consistent quantum- 
mechanical calculation either within Density functional 
Theory (DFT) or other first principles approaches, gets 
mapped into a much smaller space spanned by a set of 
atomic or atomic-like (i.e. Wannier functions) orbitals. 

The representation of the electronic structure of the ma¬ 
terials on a minimal TB basis set has been obtained with 
two main approaches: (i) The “downfolding” of the ab- 
initio electronic structure (solved in the large basis) into 
a model containing only a few bands of interest which 
are disentangled from the rest; {ii) the explicit calculation 
of the matrix elements Hap = {(j)a\H\(l>i3) using predeter¬ 
mined and fixed localized functions, typically AOs. 

In the first approach, one proceeds by selecting a sub¬ 
space B (spanned by N Bloch states |^/>„) of interest) 


of the AT-dimensional space of the solutions of the orig¬ 
inal quantum mechanical problem. The latter is found 
by representing and diagonalizing the Hamiltonian H of 
the system using a very high-quality basis set of size K, 
e.g. plane-waves with a large cut-off, a dense spatial grid, 
a large number of atomic-orbital-like Gaussian functions, 
etc. The subspace B is then projected onto a space A gen¬ 
erated by the atomic -like orbital functions, {(fa) where 
a = 1,...,M. Typically, the dimension of the B sub¬ 
space, N is much smaller of K, while the number of atomic 
orbitals is M > N and is defined by the choice of the local¬ 
ized basis set in A. Typically, in order to obtain a faithful 
representation of the electronic properties of the system in 
the smaller basis, the basis functions need to be iteratively 
optimized, thus adding a substantial computational effort. 
Implementations of this approach include: muffin-tin or¬ 
bitals of arbitrary order (A^MTO)P, maximally-localized 
Wannier functions (MLWF)P, quasi-minimal basis orbitals 
(QUAMBO)P, etc. While the optimized functions can be 
used to compute the TB matrix elements, their primary 
advantage is exploiting the information they contain to 
study the physics of the handpicked bands. For example, 
they can be used in mapping correlated bands into Hub¬ 
bard models.li^ Implementations of the second appro acltSl 
have used the non-self-consistent Harris-Foulke^^^*^ func¬ 
tional for H with the input charge density taken from the 
converged large-basis ab-initio calculation. The compu¬ 
tational bottleneck in this approach is the calculation of 
multi-center integrals (three-center and up) for the func¬ 
tional. This approach can be readily extended to find the 
charge density self-consistently, thus, allowing efficien t im¬ 
plementations of order-A^ ab-initio DFT codes.^i^^i^ 


In this work we follow the principles outlined above but 
without resorting to an explicit basis set optimization. Us¬ 
ing the eigenstates \ijjn) of H with n = 1,..., A^, one can 
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always write 

N 

Ha(3 = {(l>a\H\(j)p) « '^{(l)a\'>Pn)En{tpn\ft>p) ■ 

n—1 

Here the « sign is introduced because we restrict the sum 
to N elements (the subspace B) instead of K (the “com¬ 
plete” basis within the limits of convergence). Defining 
the matrix of overlaps Ban = {4>a\'4’n) and the diagonal 
matrix E = diag(£li, i? 2 , • ■ •, the TB Hamiltonian 

matrix, H, is expressed as 

H = BEB\ (1) 

where H is a rectangular M x N matrix. In this way, the 
computation of the TB matrix reduces to a straightfor¬ 
ward matrix operation which does not require any special 
iterative procedure as needed by some of the methods dis¬ 
cussed above in (i). This construction takes advantage 
of the full knowledge of the eigenenergies and eigenfunc¬ 
tions obtained in the large-basis calculation, in contrast 
to using the charge density only as in (ii). This scheme, 
also known as direct projection, has been tried in the past 
but was considered unreliable: even though it yielded an 
overall resemblance to the large-basis band structure, it in¬ 
troduces spurious states “randomly” scattered across the 
energy spectrum and unphysical hybridizations. See for 
instance the band structures in Fig. as well as Ref. [ 5 ] 
(Figs. 5 and 7, see also R.ef. lTBp land Figs. 1,2,4 in Ref.fT^ 
We have previously shown that accurate TB Hamil¬ 
tonians can be straightforwardly obtained from Eq. if 
only Bloch states \tpn) that project well on the selected 
AO basis set (high projectability > 95%) are included in 
the subspace B, i.e., filtering. This process introduces a 
null space, which is shifted outside the energy window of 
interest.^ In this work we present a generalized scheme 
for the construction of TB Hamiltonians in a minimal ba¬ 
sis set, suitable for cases when states with moderately 
high projectability (> 85%) are needed to be included. 
The new scheme further enhances the accuracy of the 
TB Hamiltonian and has the added advantage of making 
the TB eigenvalues insensitive to the shifting operation. 
Furthermore, we use perturbation theory to analytically 
demonstrate that the spurious states and unphysical hy¬ 
bridizations previously observed in the direct projection 
scheme are due to the presence of low projectability states. 

II. METHODOLOGY 


the projector operator P = \4‘a){4>a\ that is hermitian 

and idempotent. 

Let us consider a set of column vectors {\Bn)} obtained 
from direct projection of each Bloch wave |'0„) of the B 
subspace onto the chosen orthonormal atomic orbitals, 
P\'^n) = The elements of the vector |H„) 

are the projection coefficients B^n = The ex¬ 

plicit expression for the computation of these coefficients 
is given in Eq. |A8[ 

The projector P is represented as Pb = B^B, a. N x N 
matrix with entries 

(H'l'Hjnm = (V'nlHIV'm)- 

The diagonal elements are the “projectabilities”, defined 
as 

Pn = {B^B)nn = {Bn\Bn) = {tpn\P\i’n}, 

that measure to what extend the eigenstate ijn is well- 
described in the space A specified by the projector P. 

While by construction tr(P) = M, the trace of the ma¬ 
trix Pb in B becomes 

N 

ir{PB) = Y,{'4’n\PHn) < N < M. 

n—1 

If A is complete, all i/'n states project perfectly (p„ = I) 
and the trace equals N, the size of B. The deviation from 
Y is a criteria to assess the accuracy of the TB represen¬ 
tation in A with respect to the electronic structure in the 
B subspace. 

High projectabilities are expected for the states in the 
lowest bands and poor projectabilities at higher energies. 
One would expect to obtain accurate TB eigenvalues and 
eigenvectors for the states with the largest « I; how¬ 
ever, as proved in Sec. |TH] states with low projectability, 
when folded in the TB Hamiltonian, hinder the accuracy 
of the results. 

In order to represent well the electronic structure of the 
system one needs to exclude from B the “bad” states with 
low projectability by choosing N accordingly (p„ larger 
than a chosen threshold for each n = 1,...,Y). The 
procedure involves the construction of the tight-binding 
Hamiltonian following Eq. with a normalized set of col¬ 
umn vectors |A„) = \Bn )/such that (A„|A„) = 1. 

These states are used to build the initial TB Hamilto¬ 
nian: 


In our work we use plane-waves (PW) as the large ba¬ 
sis for the ab-initio calculation of the Bloch states 
The wave vector k index is suppressed so the analysis for 
a periodic system can be understood to be at a particu¬ 
lar k point. We choose dim(A) = Matomic-like localized 
orbitals \(j>a) (with M > N = dim(H)) which we assume 
to be an orthonormal set {4>a\4>i3) = ^a/S- These could 
be Wannier functions or, more pragmatically, Lowdin or¬ 
bitals. The restriction to this subspace is obtained through 


H = AEA ^, 


( 2 ) 


where the N columns of the matrix A are the vectors |A„). 
This product is an M x M matrix constructed using only 
N states. Because of this construction, H is singular with 
an unphysical null space, Af, of size M — N that com¬ 
promises the accuracy of the eigenvalues (see Sec. HI). In 


order to remove the effect of the null space we perform an 
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orthogonal projection (see Ref. Uni using the set of vectors 

{|^n)} 


Q^ = Im-A{A^A)-^A\ (3) 

where Im is the M x M identity matrix. Reconstructing 
the TB Hamiltonian as 

H^ = H + (4) 

is possible to shift the eigenenergies corresponding to the 
null space elements to an arbitrary energy k, away from 
the band with good projectability. In practice, if only very 
high-projectability states are considered, A^A is close to 
the identity, and the shifting matrix can be approximated 
by Qm ~ Im — AA\ avoiding the matrix inversion. This 
approximation introduces a small k dependence into the 
states of the TB subspace. This dependence can be safely 
neglected when using a very high-projectability filtering 
criteria {e.g. Pn > 0.95 in Ref. mi) or applying only small 
values of k, otherwise the exact expression in Eq. is 
required for a faithful description of the energy bands. 


III. THE EFFECT OF LOW PROJECTABILITY 
STATES IN THE TB HAMILTONIAN 


Minimal basis set have proved satisfactory to achieve 
accurate TB matrices for periodic systems using the fil¬ 
tering procedure. However, if more unoccupied bands of 
high projectability are needed for a particular application, 
one can achieve that by progressively increasing the size 
of the AO basis set, e.g. from single zeta (SZ) (minimal) 
to double zeta (DZ), etc. effectively increasing the size of 
Nonetheless, for most cases it is more advantageous 
to trade-off some accuracy away from the Fermi energy 
for the convenience of still dealing with TB matrices of 
smaller sizes, i.e. to keep the basis minimal, especially in 
the study of systems with large number of atoms. This 
can be achieved by including bands with moderately high 
projectability (p„ > 0.85) in the construction of the TB 
matrix. However, as we discussed in Sec. |llj bands with 
low projectability affect the accuracy of the TB represen¬ 
tation in A. 

To learn about the eigenvalues of we start by ap¬ 
plying it to A. One gets 

H^A = HA + nQjg-A = HA = AEP , 

where P = A^A. To find an analytical expression for 
the eigenvalues of we assume the number of states 
ifin to be equal to the number of AOs {N = M), so that 
A is square and invertible.^ Then, the expression above, 
Hk,A = A{EP), has a M x M square matrix EP which 
is the representation of H in some linearly independent 
basis (columns of A) and its eigenvalues are also those of 
Hk (and H) that we call E. 


EP, with the diagonal pulled out as a perturbation is 
then: 


EP — diag(Pii£'i,..., PmmEm)+ 

I 0 P12E1 P13E1 

P21E2 0 P23E2 

P31E3 P32E3 0 


\ 


: Pm-i,mEm-i 

\ Pm,m-iEm 0 j 

First of all, if any = 1, then that column and row of 
the perturbation are zero. The diagonal element is decou¬ 
pled from the rest of the matrix so one eigenvalue will be 
exactly = E^. In other words, a perfect representation 
of the exact wavefunction means the basis is complete for 
that state and that the action of H in that basis will be 
perfectly described. 

Second, if one ignores the off diagonal entries, it can be 
seen that the eigenvalues of would be scaled versions 
of H where each eigenvalue is being scaled by its pro¬ 
jectability so the energies of Hk would be « PnnEn. 
Therefore, bad projectability will incorrectly deliver a TB 
eigenvalue close to zero. This has been the source of much 
trouble in previous TB methods without filtering. 

Third, the off-diagonal elements of the perturbation ma¬ 
trix lead to hybridization. Namely, the lack of perfect 
projectability leads to level repulsion and further changes 
of the TB eigenvalues, beyond the scaling by men¬ 
tioned above. The perturbation matrix is small, since (i) 
the off-diagonal elements of P are negligible for all states 
with high projectability; and (ii) they are still consider¬ 
ably small for the states with lower projectability (See 
Appendix Second order perturbation theory on the 
off-diagonal elements gives the analytical expression for 
the n-th eigenvalue of H, which is called E^, changing 
from PnnEn tO 


En — Pn.7i.En 


E 

3¥=n 


\P,n\'^EnEj 

PnnEn PjjEj 


0{P^ 


( 5 ) 


In addition to showing that the hybridization changes the 
energies, the formula shows that the changes due to hy¬ 
bridization can be much larger than one would naively ex¬ 
pect based only on looking at the small Pjn entries: first 
the numerator has an additional Ej energy factor, and sec¬ 
ond, the energy difference in the denominator is based on 
the scaled energies which means the energy difference can 
be smaller than between the actual eigenvalues (especially 
when both and Pjj are significantly smaller than one) 
thus enhancing the contribution of the hybridization. 


IV. RESULTS AND DISCUSSION 

Here we illustrate the effects of low-projectability bands 
on the accuracy of the TB matrices derived for a molec¬ 
ular system, benzene, and a crystal, cobalt antimonide. 
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The ab-initio calculations were performed with plane-wave 
DFT codes: VASr^^ and/or quantum espressoP^ 

A. Benzene 

In the case of benzene we computed the electronic 
structure using VASP within the projector-augmented- 
wave (PAW) methocP^ and the Perdew-Burke-Ernzerhof 
(PBEjpSl functional. We determined the molecular wave- 
functions (molecular orbitals) of an isolated benzene 
molecule in a large cubic supercell of 15 A side using P- 
point to sample the reciprocal space and a kinetic-energy 
cutoff of 29.4 Ry. 

The detailed procedure to compute the projection coef¬ 
ficients, Ban, is discussed in Appendix [A| The states 
are projected onto a minimal basis set of M = 30 AOs 
(C:2s, 2p; and H:ls) taken from public repositories.^ 


11 17 22 



FIG. 1. (Color online) Projector matrix |Pb| of benzene on a 
minimal AO basis set for the 23 molecular orbitals of lowest 
energy. The diagonal elements are the projectability numbers 
p„. The presence of non-zero off-diagonal elements {Bm\Bn} 
(m yf n) reflects the non-orthonormality of the vectors |Bn). 


The diagonal elements of the projector matrix Pg = 
B shown Fig. are the projectability numbers for 
each Bloch state t/jn- In the chosen AO basis set, the 
17 Bloch states of lowest energy have high projectability 
{pn > 0.88) whereas states 18 < n < 20 have low pro¬ 
jectability (0.20 < Pn < 0.32). Moreover, higher-energy 
states (21 < n < 23) are not projectable in this particu¬ 
lar AO set {pn = 0.1,0.005,0.076, respectively.) As dis¬ 
cussed above, an accurate TB Hamiltonian matrix can be 
constructed by filtering out states with low projectabil¬ 
ity. Therefore, considering only the lowest iV = 17 states 
{pn > 0.88) yields TB eigenenergies that are in excellent 
agreement with the DFT values. The black dots in Fig. 
for N=17 show a maximum deviation from the DFT en¬ 
ergies (gray lines) of only 5 meV. 

In order to study the effect of states with lower pro- 
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FIG. 2. (Golor online) (a) Evolution of the tight-binding 
eigenvalues with increasing number of low-projectability Bloch 
states. N is the number of states used in the construction of 
the TB matrix, where the first 17 states are those of high pro¬ 
jectability. In all case, the null states (not seen) are shifted by 
K = 8 eV. (b) Zoom-in around the low-projectability TB eigen¬ 
values. The red dots mark the TB eigenenergy corresponding 
to the 18th Bloch state .(c) Energy variation of each TB state 
n against increasing N. The reference energy E/ corresponds 
to the Hamiltonian without low-projectability states (A=17). 


jectability on the accuracy of the TB eigenenergies, we 
intentionally relax the filtering criterion to progressively 
include some states with lower projectability {N > 17) in 
the B subspace. The vectors |A„) corresponding to the 
low-projectability states are left unnormalized, i.e., 

lAX_j\Bn)/y/p^, if Pn>0.85 
" [ |E„) otherwise. 

The normalization of the high-projectability vectors |A„) 
artificially makes the corresponding diagonal elements 
{PB)nn = Pn (~ 1) equal to 1, but this small change does 
not alter the analysis that follows. 

Expectedly, the TB eigenenergies En corresponding to 
low-projectability states largely underestimate the DFT 
values En- Those TB energies are seen in the zoom-in box 
in Figs. and b. 

The inclusion of the lower projectability state n = 18, 
for instance, yields the “bad” TB eigenvalue Ais = 0.9732 
eV. The large underestimation with respect to Eis, seen in 
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Fig.§ is accounted by directly scaling the DFT values by 
the projectability, i.e. Pigps-E'is = 0.9903 eV, as discussed 
before. The evolution of Sig with increasing N is shown 
in Fig. using red dots for visual aid. 

Including the degenerate states n =19,20 also yields 
strongly underestimated values Eig = E 20 = 0.7621 eV 
and in agreement with Pig^igEig = 0.7877 eV. 

Moreover, including the states 21 < n < 23, which 
have even smaller projectabilities, yields TB eigenvalues 
close to zero: 0.0, 0.0, 0.2949 eV, respectively. Ex¬ 
cept for E 21 , these values compare well to = 

0.3984,0.0218,0.3117 eV, respectively. The departure of 
E 21 from P 2 i, 2 iE 2 i is due to hybridization effects discussed 
later in the text. 

As hinted by Eq. a state |A„) hybridizes with an¬ 
other states \Aj) via a non-zero off-diagonal element of 
Pnj = {An\Aj). The two diagonal blocks, of size 17 x 17 
and 6x6, seen in Fig. correspond to high and low pro¬ 
jectability states. The 17 x 6 and 6x17 off-diagonal blocks 
allow hybridizations between both types of states. Since 
all elements in the off-diagonal blocks are small (< 0.0729 
eV), the low-projectability states are expected to have only 
a small impact on the high projectability ones. This is 
confirmed by inspecting the variations of the “good” TB 
eigenenergies (E„, n = 1-17) while increasing the size N 
of the subspace B. The variations are small and not notice¬ 
able in Fig. Instead, we plot the TB energies relative 
to reference values A* in Fig. [^. E* are the TB eigenen¬ 
ergies when A^=17, that is, the high-projectability case. 
The maximum variation found is 21.8 meV and happens 
for Ell in the fifth panel (once state n=22 is included). 
This is consistent with the maximum element of the off- 
diagonal block happening at |T’ii, 22 | = 0.0714. We find 
that Pi 1,22 is the only non-zero off-diagonal element in 
the 11th row (and column) of the projector matrix, the 
energy variation can be directly attributed to the overlap 
between |Aii) and IA 22 ) following the hybridization mech¬ 
anism depicted by the yellow circles and arrows in Fig. [l] 
The variation is well estimated by the second-order per¬ 
turbation model, Eq. which reduces to: 


APii 


|Pll,22pPllp22 
PliqiPll ~ P22,22p22 


= 21.7 meV. 


All hybridizations due to the off-diagonal elements in 
Fig. 0 (or similarly, of the matrix P) translate into peaks 
in Fig. The number of peaks in each panel increases as 
more low-projectability states are progressively included. 
Every new peak n that appears in a particular panel 
N reflects the hybridization between a low-projectability 
state—the one newly introduced in panel N —and the 
“good” TB eigenstate n. Each new peak can be directly 
traced to a non-zero off-diagonal element in Fig. For 
instance, the peaks at n=l,6 in panel N=18 are due to 
Pi, 18 , Pe.is; peaks n=3,10 in panel 77=19 to Ps.ig, Pio.ig; 
etc. 


The elements of the off-diagonal blocks yield only small 
fluctuations; however, the overall maximum off-diagonal 


element |Pi 8 . 2 i| = 0.1757 is inside the smaller 6 x 6 di¬ 
agonal block. This indicates that hybridizations between 
the low-projectability states jAig) and IA 21 ) are stronger. 
Hybridization causes the level repulsion of Pig (red dot) 
along the upward arrow observed in Fig. [2> at the in¬ 
troduction of n=21. The repulsion shifts up the level by 
0.3953 eV from Pig^igPig, which is the value expected in 
the absence of hybridization, marked by the red cross. The 
leading hybridization mechanism is depicted using green 
circles and lines in Fig. The perturbative estimate of 
the level repulsion is 


APig 


|Pl8,21pPl8p21 

P18,18^18 ~ P21,21p21 


= 0.6656 eV. 


P 21 has the opposite level repulsion AP 21 = —APig 
(downward arrow), which explains the discrepancy 
between the value in the absence of hybridization 
P 21 . 21 P 21 = 0.3984 eV (black cross) and the actual TB 
energy P 21 « 0 eV discussed earlier. 
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FIG. 3. (Color online) Behavior of the null and low- 
projectability TB eigenstates under the shifting operation using 
(a) the exact Qjv = Im — A{A^A )~^or (b) the approxi¬ 
mated Qj\/ ~ Im — AA^ shifting matrix. 


The construction in Eq. introduces a null space A/” 
containing {M — N) degenerate eigenenergies Ej\f = 0 eV. 
The matrix Qj^ = Im — A{A^A )~^is used to selectively 
move the null subspace upwards in energy by the control 
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parameter k, without affecting the remaining TB energies. 
With the exact Q^f the TB values do not acquire a depen¬ 
dence on K. The evolution of the TB eigenvalues with 
K in Fig. ii readily shows that only the degenerate null 
eigenenergies have a dependence on k (marked with solid 
black lines). The low-projectability TB states (dots about 
0.85 eV in panels A^=18 and 7V=19) do not belong to the 
null space and therefore are also independent of k. As ar¬ 
gued in Sec.im the shifting matrix can be approximated by 
Qm ~ Im — AAl, which avoids a matrix inversion, but at 
the cost of introducing a small k dependence to the good 
TB values. The approximation is safe when using a high- 
projectability filtering criterion {N=17). In this case both 
the exact (first panel in Fig. [^a)) and the approximated 
Qj^ (first panel in Fig. |^b)) yield the same TB eigenval¬ 
ues. Nonetheless, the energy deviation due to the k de¬ 
pendence introduced by the approximated Qj\f can become 
significant when low-projectability Bloch states are intro¬ 
duced, for instance, compare the second (and third) panels 
in Figs. da) andj^b). With the exact Q^/ (Fig. [^a)) the 
low-projectability TB eigenvalues around 0.85 eV remain 
flat whereas they acquire a chiefly linear k dependence 
(marked with blue lines) when using the approximated 
(Fig.db)). 


(a) (b) 


44 55 58 



0.00 0.25 0.50 0.75 1.00 


30 35 40 45 50 55 
Bloch state n 

FIG. 4. (Color online) (a) Minimum projectability per band 
n over all k points for CoSbs on the chosen minimal basis set. A 
sharp decline of the projectability {p„ < 0.2) is seen for states 
n > 58. (b) Projector matrix \Pb\ at ko = (0.4,—0.4, 0.5), in 
reciprocal coordinates. Only the matrix elements from 44 < 
n < 58 are shown. 


B. CoSbg 

We analyze CoSba as an example of a periodic solid. 
CoSba is a typical binary skutterudite compound with cu¬ 
bic structure and space group #204. Skutterudites are 
among the most promising thermoelectric materials.^ 


We use the quantum espresso suite of ab-initio codes 
to obtain the Bloch wavefunctions and the matrix B of 
projection coefficients. The wavefunctions are obtained 
using the PBE functional with an energy cutoff of 50 Ry 
and the PAW data set from the PSlibrary 1.0.0.^ We 
choose a minimal basis set to project onto composed of 
M = 84 AOs (Co: 4s,4p, 3d and Sb: 5s, 5p, taken from 
the PAW data set). 

The projectability of all occupied DFT bands (n < 48) 
is very high > 0.97 and progressively decreases for 
the unoccupied bands at higher energies (see Fig. |^a)). 
The basis set supports 9 unoccupied bands with high pro¬ 
jectability of pn > 0.87 (49 < n < 57) before declining 
to poor values of < 0.15 for n > 58. An accurate 
TB Hamiltonian is obtained when considering only Bloch 
states of high-projectability (N=57} as confirmed by the 
excellent match between the TB (black and red dots) and 
the DFT bands (gray lines) seen in Fig. [^a). The band 
n=55 is shown in red. 


To illustrate the effect of hybridization on the accuracy 
of the TB Hamiltonian, we include one Bloch state n=58 
of low projectability (psg = 0.15) in the construction of the 
Hamiltonian. The corresponding TB band (not present in 
panel a) is shown in Fig. in blue. 

First, it is seen that the TB band n=58 (blue dots in 
panel b) does not reproduce the reference DFT band, es¬ 
pecially in regions of reciprocal space with the lowest pro¬ 
jectability (around H) where the TB eigenvalues strongly 
underestimate the DFT values and are consistent with be¬ 
ing scaled by their projectabilities, i.e. Pnn^n ~ 0.4 eV. 

Second, the hdelity of the band n=55 (red dots) is no¬ 
ticeably reduced around H, with respect to panel (a), due 
to hybridization with the low-projectability band n=58 
(blue dots), which leads to level repulsions. The in¬ 
set in Fig. shows the repulsion of the eigenstates 
at ko = (0.4,—0.4, 0.5), in reciprocal coordinates. The 
crosses mark the values of #55 (red) and Ess (blue) ex¬ 
pected in the absence of hybridization (« PnnEn). The 
level repulsion due to hybridization is shown along the 
arrows. As seen in Fig. [^, the predominant non-zero off- 
diagonal element at ko is ^ 55 , 58 , which indicates that the 
hybridization primarily involves only IA 55 ) and [Ass) as 
depicted by the green circles and lines in Fig. There¬ 
fore, the repulsion of # 55 —the magnitude of which is given 
by the size of the arrows in the inset—can be analytically 
estimated by the perturbation formula as 


AEss = —AEss 


\P55,58\‘^E55E58 

Pbb,bbEsb — Pb8,b8Es8 


0.1928 eV, 


which in agreement to the actual value of 0.1381 eV. 


All the eigenenergies seen in both panels (a) and (b) of 
Fig.[5l *.e. the TB subspace, are insensitive to any cho¬ 
sen value of K since the exact shifting matrix Qj^f is used. 
The null subspace (eigenvalues not seen) have been rigidly 
shifted by K = 2 eV outside the region of interest. In panel 
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(b) (c) 



FIG. 5. (Color online) Tight-binding eigenenergies for CoSbs. (a) Accurate TB bands are obtained when adopting a high- 
projectability filtering criterion (p„ > 0.87) in the TB Hamiltonian. The band n=55 is shown in red. (b) The inclusion of a 
low-projectability band (psg = 0.1474) introduces the eigenenergies seen in bine, which induces unphysical hybridizations with 
the band n=55. The arrows in the inset, whose size is lAJSssI, illustrates the level repulsion at ko = (0.4, —0.4,0.5), in reciprocal 
coordinates, due to hybridization, (c) The TB Hamiltonian is built with the approximated shifting matrix Qj^ and k = 1.0 eV 
to show the distinct dependences on k of the eigenstates. The reference DFT bands are shown with gray lines. Brillouin zone 
integration follows the AFLOW standard as discussed in Ref. m 


(c) we re-compute the TB Hamiltonian from (b) but us¬ 
ing the approximated shifting matrix {Qj^ « Im — AA^) 
and K = 1.0 eV. Different shifting patterns are observed: 
(i) The eigenvalues of the null subspace shift rigidly with 
the value of k and, thus, are pinned along the horizontal 
line at 1.0 eV. {ii) The “unhybridized” bands of the TB 
subspace (black dots) show no noticeable difference with 
respect to panel Fig. Ep- This confirms that while they for¬ 
mally acquire a n dependence, introduced by the approx¬ 
imated Qf^f, the effect is negligible for bands with high 
projectability. (m) The high-projectability band n=55 
(red dots) should also be insensitive to k; nonetheless, it 
acquires a more noticeable dependence indirectly via its 
hybridization to the K-dependent low-projectability band 
n=58. Consequently, the most noticeable changes of this 
band with respect to (b)P^ happen around H where the hy¬ 
bridization is stronger, (iv) The low-projectability band 
n=58 (blue dots) shows a noticeable dependence on k, 
especially around the lowest-projectability k points. For 
instance, the value of k = 1.0 eV effectively shifts the 
states around H by ^ 0.9 eV, i.e. from ^ 0.4 eV [as in 
Fig. [^b)l^ to ~ 1.3 eV. 

V. CONCLUSIONS 

In this paper we have outlined a noniterative scheme to 
derive highly accurate ab-initio TB Hamiltonian matrices 
in a minimal basis set representation. 

Minimal basis sets may be insufficient to converge self- 


consistent quantum-mechanical calculations with linear 
combination of atomic orbitals, however, they are ade¬ 
quate for the purpose of projecting wavefunctions obtained 
with fully converged basis and building the reduced TB 
matrices. 

Low-projectability Bloch states have spurious effects 
when included in the construction of the TB matrix. We 
have unambiguously shown the underlying mechanism for 
the formation of the spurious states. The removal of those 
states, via the application of a shifting matrix, delivers 
accurate TB matrices. 

We have introduced an expression for the shifting ma¬ 
trix of a nonorthogonal set of vectors. This expression 
improves the quality of the TB Hamiltonian by removing 
any unwanted dependence that the shifting procedure had 
on the TB eigenstates of interest. 

Appendix A: Projection of Bloch states on pseudo 
atomic orbitals 

The pseudo-wavefunction Bloch state \4’nk) is expanded 
in plane-wave basis jk -|- G) as 

|^„k) =^GG„k|k + G). (Al) 

G 

The plane-wave basis 

(rjk + G) = (A2) 

V ^0 

















is defined to be normalized to 1 over the volume of the 
primitive unit cell rio- The orthonormality of the basis 
(k + G|k + G') = Jgg' allows the expansion coefficients 
to be defined by the projection 

C'cnk = (k + G|'!/j„k) (A3) 

In the US/PAW pseudopotential formalisms, the pro¬ 
jection of the all-electron (AE) wavefunctions V'nk onto 
an atomic orbital (f)^ is computed in terms of their corre¬ 
sponding pseudized quantities ^„k and and the overlap 
operator ,§ = i -|- \l3\)Qlj{l3'}j\- The pseudo atomic 

orbitals (PAO) and beta projectors are defined as 

U^) = Rl{r)Yr{v) (A4) 

where /i = {Ilm\ is a composite index of the ion center I 
and quantum numbers {lm\ of the PAO. The real-space 
beta projectors are analogously definecP^, with i = {Im}, 
as: 

Pu(v) = RUr)Yr{r) (AS) 

The localized basis |0j^) for periodic calculations is con¬ 
structed from Bloch sums of the PAOs. 

R 

where N is the number of lattice vectors R. The Bloch 
sum for the AO basis |(/>J^) follows the same definition. 
Notice that the factor implies normalization of \(j>^) to 
1 over the primitive unit cell, which is consistent with the 
normalization of the plane-wave basis in Eq. |A2| 

Then, the projection coefficients are calculated in term 
of the pseudized quantities: 

= (<('}) I V'nk) = 

= (<^;)|^nk) + E ■ (A7) 

GG' 

lij 

The integrals are more efficiently computed in the 
|G) basis. Using the identity i = Vg |GI')(G| = 
X:g |k-f G)(k-HG| one has 

E(<^)^|k + G)(k-bG|^„k) + 

GG^ 

^ (<^;:|k + G) (k + G\Pl)Ql^ {Pi |k -f G')(k + G>„k) 

GG' 

Hj 


where the objects in brackets are given in Eqs. B6 and B7 


The coefficients expand the Bloch state on a PAO 
basis {(j>fj,}- Furthermore, the coefficients on a Lowdin 
orthonormal basis are readily obtained by 


= (A8) 

i 

where the upper bar symbol is used to indicate orthonor¬ 
mality and = {PIPI) is the matrix of overlaps be¬ 
tween PAOs. 

Lowdin coefficients and orbitals are assumed through¬ 
out the main text where we drop the upper bar and k 
superscript in the notation of B^. 


Appendix B: Projection of the plane-waves basis on 
pseudo atomic orbitals 


Using the relation 1 = / |r)(r|dr and Eqs. 


A2 


and 


A6 


to evaluate the projection (0|^|k - 


G) one has 


(^;)|k + G) 


dr(0{^|r)(r|k-l- G) 


(Bl) 


iV-^/Og 




R 


(B2) 


oi(k+G).Tj, 

N y/np) 


„i(k+G)- 


R 

(B3) 

/dre*('"+®)"^7(r). (B4) 


Using the plane-wave expansion gdk+G).r _ 

Ti'm' 4^f''iz'(|k + G|r)r-'*(irrG)y;U'(r), where 

the hat notation in dica tes the directional angles of the 
vector under it; Eq. |A4[ and dr = smOdrdOdip, the last 
expression reduces to 


-1 4^gd^t-G)-r^ ^^ - 

(^l^lk + G) =-*(k + G) 

X J r^R’ll{r)ji'{\k + G\r)dr j Y,™*(r)r,U'(r) sinddddy:. 

(B5) 


Using the definition in Eq. |A3| for the expansion coeffi¬ 
cients: 

Bln = {4>l\^n^.) = + G)CG„k + 

G 

^ |k -y G) (k + G\Pl)Ql iPl |k -f G')CG-„k , 

GG' 

Hj 


With the the normalization identity 
J Y/^*{9, (l))YiY {0,(j)) sin 9dOdip = SwSmm', we arrive 
to the final expression: 

(^;j|k + G) = /G^kFr*(l^) J r^Rtir)ji{\k + G\T)dr 

(B6) 


with /Gpk = dTTi'Uo " 
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Analogously for the projection on the beta functions: 

(/3]^Jk + G) = f r^R^j,ir)ji{\k + G\r)dr . 

(B7) 

Appendix C: OfF-diagonal elements of the projector 
matrix 

Given the projector matrix 

Prim = {B^B)nm = {'^n\P\i’m) , 

since ipn forms a complete Hilbert space, the matrix P will 
also be a projection operator by closure. Namely, 

Prim — (,P ')nm — ^ ^ PnjPjm j 
j 

so there is a constraint for the important diagonal elements 
(and using the Hermitian nature of the matrix P) : 

Pn — Pn P 'y ^ I Prim \ ■ 

m^n 

This expression puts an upper bound of 1 on the diago¬ 
nal elements. One can also define an upper bound on any 


off diagonal element via 

\Pnm\ < min(Vp„ -pI, Vpm - Pm) ■ 

Notice that if = 1 (or « 1) which means perfect 
projection, then Pnm = Pmn = 0 (or « 0) Vm ^ n so the 
entire n-th column and row of P is zero (excluding the 
diagonal which is 1). 

For cases of smaller projectability < 1, each off di¬ 
agonal entry will still be much smaller than since there 
are many of them in the sum rule; but the sum of their 
squares must add up to — Pn- 
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